arXiv:1508.02644v2 [cond-mat.stat-mech] 16 Dec 2015 


Crossover from ff to a relaxation in cooperative facilitation dynamics 

Mauro Sellitto 

Dipartimento di Ingegneria Industrials e dell’Informazione, Seconda Universita di Napoli, 

Real Casa dell’Annunziata, Via Roma 29, 1-81031 Aversa (CE), Italy 

/? and a relaxation processes are dynamical scaling regimes of glassy systems occurring on two 
separate time scales which both diverge as the glass state is approached. We study here the crossover 
scaling from (3- to a- relaxation in the cooperative facilitation scenario (CFS) and show that it is 
quantitatively described, with no adjustable parameter, by the leading order asymptotic formulas 
for scaling predicted by the mode-coupling theory (MCT). These results establish (i) the mutual 
universality of the MCT and CFS, and (ii) the existence of a purely dynamic realization of MCT 
which is distinct from the well-established random first order transition scenario for disordered 
systems. Some implications of the emerging kinetic-static duality are discussed. 


The glass state is obtained with striking simplicity in 
a variety of substances pQ. Yet, its fundamental nature 
is one of the most enduring puzzles of condensed matter 
physics. From a dynamical point of view, what makes the 
glass relaxation so peculiar is the existence of two time 
scales which both grow dramatically as the separation a 
from the glass state vanishes. The first scale, t a , refers 
to the rattling motion of particles in the cage formed by 
their neighbors (/3 relaxation). The second one, t' a t a , 
is generally ascribed to the slower cooperative rearrange¬ 
ments of cages that allow long-term particle diffusion (a 
relaxation). This two-step relaxation pattern was first 
predicted by mode-coupling theory (MCT) (2], and en¬ 
tails as precursors phenomena strongly temperature or 
density dependent spectra exhibiting relaxation stretch¬ 
ing. Pioneer work first identified these glass precursors in 
colloidal suspensions of hard spheres [3] and in molecular- 
dynamics simulations for a Lennard-Jones mixture |4j. 

The crossover from ft to a relaxation, which occurs for 
times large on scale t a but small on scale tfi, constitutes a 
crucial test for our understanding of glassy dynamics i5]. 
It is remarkable because the ratio t' a /t a is itself diverging 
as a —> 0. At the same time, however, its actual investi¬ 
gation is delicate because of the sensitive dependence on 
three quantities: the separation parameter a, the system- 
specific exponent parameter A, and the arrested part of 
correlation, /. In realistic systems these quantities can 
be inferred only with limited accuracy (or may even be 
not well defined). This makes it difficult to assess unam¬ 
biguously the MCT predictions and only a few cases have 
been analyzed 00. In this Letter we provide a test of 
the crossover scaling in the cooperative facilitation sce¬ 
nario (CFS) 0. Several results have already suggested 
a close analogy between MCT and CFS Here we 

show that this relation is, in fact, quantitative and deep. 

MCT scaling laws. - Without loss of generality we 
refer to the ideal version of MCT, in which the normalized 
density-fluctuation correlator $(t) has no wave-vector 
dependence (for an exhaustive description of MCT, see 
Ref. j2j). Near the glass bifurcation singularity the cor¬ 
relator exhibits a plateau at some f c . Near this plateau 
the small-er dynamics on scale t a is ruled by the first 


scaling law: 

$ 0 ) - fc = c a g(i), i=l~. ( 1 ) 

Here the correlation scale c a and the time scale t„ read 
c <7 cx y/ft, t& oc |cr|~«, (2) 


where a is some critical exponent. The a -independent 
master function g(t) obeys the scale-invariant equation 


- 1 + A git) 2 - d f f git - t') g(t') d t' = 0, (3) 

Jo 

to be solved for the initial condition lining g(t)t a = 1. 
The exponent parameter, A, quantifies all properties of 
git). In particular, it determines the exponent a via 


r(i-q ) 2 

T(1 - 2a) 


(4) 


Two time regions can be distinguished in the first scaling 
regime where, depending on whether $ is above or below 
f c , the master function g{t) can be approximated as 

f gait) = t~ a - A\ t a , for <F > f c 

g(i)^{ „ (5) 

{ Obit ) = ~B t b + ^ t b , for $ < f c 

with an error decreasing faster than t a in the former case 
and t~ b in the latter. The coefficients A\ and B\ read 

^-=T(l-a)r(l + a)-A ( 6 ) 

±-=Til-b)Til + b)-X, ( 7 ) 

while B is a positive quantity fixed only by A. The domi¬ 
nant contribution to gbit), —Bt b , is called von Schweidler 
decay. The von Schweidler exponent b is determined by 


T(1 + b) 2 

r(l + 26) 


( 8 ) 


with 0 < b < 1. The function git) exhibits a zero at 
some t*. The part for t preceding t* deals with deviation 
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FIG. 1. Correlation decay for temperature T £ [0.49, 10] FIG. 2. MCT first scaling law in the CFS. Rescaled correlator 
in the CFS on a Bethe lattice with coordination 2 = 4 and vs reS caled time variable 't/C* =t/t*. 
facilitation parameter f = 2. System size N = 2 23-24 . 


from the critical decay. The part for t > t* describes the 
approach toward the von Schweidler decay. The quanti¬ 
ties t* and B have to be calculated by solving Eq. ^ . 
The numerical solution [5] shows that there is an interval 
for t centered by t* , where the three functions g a , gb and 
g are very close to each other. Therefore, to locate ap¬ 
proximately t* and B we set g a {t*) = gt{t*) = 0, which 

gives t* ~ A 1 1 ,/2a and B ~ \JA^ a B±. These approxi¬ 
mations permit us to base the comparison with the CFS 
on elementary formulas. 

Substituting the von Schweidler decay into Eq. |l]) , one 
gets the von Schweidler decay law 

m-fc = -(p) , ( 9 ) 

where the new relevant time scale is now 

CcxM-, 7=^- (10) 

von Schweidler’s law describes the small-u dynamics for 
times intermediate between two diverging time scales, 
1, <C t « It is important because of the connec¬ 
tion it establishes between the late part of /3 relaxation 
to the early part of the a relaxation. The latter deals 
with the plateau-below dynamics on scale t' a and is glob¬ 
ally governed by the second scaling law (also known as 
the time-temperature superposition principle): 

= ( 11 ) 

Here $ is another er-independent master function. As 
an example, Figs. m and [3] show the typical two-step 
relaxation and the associated MCT scaling laws for the 
CFS that we now describe. 


CFS. - In the dynamic facilitation approach the meso¬ 
scopic structure of a liquid at temperature T is rep¬ 
resented by an assembly of binary spins, Si = ±1, 
i = 1, ..., N, whose value depends on whether the lo¬ 
cal liquid density at site i is higher or lower than the 
average. Energetic interactions are absent but the spin 
dynamics is facilitated (or constrained) |8j: At each time 
step a randomly chosen spin Si is flipped with probabil¬ 
ity w(si —> —Si) = min {l, e -Si/,T }, provided that the 
spin Si is surrounded by at least f nearby (liquidlike) up 
spins. This constraint mimics the cage effect: when the 
temperature is low enough the fraction of liquidlike spins 
is vanishingly small and, consequently, spin relaxation 
may involve a large number of cooperative spin flips over 
regions with increasing size. Here, the relevant range of 
the facilitation parameter, f, is 1 < f < z — 1, where z 
is the lattice coordination. The correlator of our interest 
shall be the persistence $(£), i.e., the probability that a 
spin has never flipped between times 0 and t. Its long¬ 
time limit, the probability that a spin is permanently 
frozen, is precisely the analog of the nonergodicity pa¬ 
rameter /. On the Bethe lattice this cooperative dynam¬ 
ics exhibits features qualitatively similar to MCT [9HI33, 
though one would, naively, expect that the conventional 
form of MCT fails for this type of model as static correla¬ 
tions vanish. In fact, previous works on related facilitated 
systems on finite dimensional lattices showed that MCT- 
like approximations are unable to describe the overall 
glassy dynamics [T4U16] . This can be generally under¬ 
stood through the connection of CFS with the bootstrap 
percolation transition which is smeared out in finite di¬ 
mensions. For this reason, our evaluation of the MCT 
status is carried out on the Bethe lattice, which is the 
first natural step of a statistical mechanics treatment. 
This is also relevant because: (i) MCT scaling regimes 
are sometimes hardly observed in numerical simulations 
of disordered systems [201 2T] , and (ii) MCT shows a be- 
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FIG. 3. MCT second scaling law in the CFS. 


havior which, in the limit of large space dimensionality, 
does not conform to the replica theory mm- 

To compare CFS with MCT we have first checked the 
/3- and a-scaling laws along with Eqs. ([2]) and (10). We 
proceed as follows. The time scale t a is estimated as 


t a = f [$(£) - $(**)] d t, (12) 

Jo 

where t* is the time to cross the plateau, $(£*) = f c . 
This is easily done as f c is known exactly. For the Bethe 
lattice we consider here, z = 4 and f = 2, we have 
T c = 0.48089... and f c = 0.67309...[5]- Figure |4| shows 
that t a obeys the power law t a oc a~ ^ with an exponent 
a — 0.27. This value is obtained by a fit over the temper¬ 
ature range [0.52, 0.82]. Notice, that a is the only critical 
exponent we estimate numerically here. The correlation 
scale c a , instead, is obtained from the MCT relation |2J: 

c a = Vl^X(f- f c ). ( 13 ) 


Here the jump of the order parameter, / — / c , is deduced 
from the exact calculation on the Bethe lattice [5] , while A 
is estimated by exploiting Eq. Q which gives A ~ 0.815. 
The result for the correlation scale c a is shown in Fig. [4] 
along with the MCT prediction. The excellent collapse of 
the rescaled relaxation data showed in Fig. [2] is obtained 
by using the c a calculated in this way. We then get b ~ 
0.45 and 7 ~ 2.96 from Eqs. (| 8 ]) and ( fl0| ). This latter 
value is consistent with that found in Ref. 0 (7 - 2.9). 
As a consistency check we then estimate the time scale t' a 
through Eq. (12) with t* such that $(£*) = 0. Figure [4] 
shows that and the ratio t' a /t a are correctly predicted 
by MCT. From t' a we finally get the second scaling law for 
the a-relaxation process that we find to hold in a wide 
range of temperatures above T c , see Fig. [3} 

We now consider in detail the MCT predictions for the 
asymptotic leading order corrections to the critical de¬ 
cay laws, Eq. ([ 5 ]). The coefficients entering Eq. |5]) as 


FIG. 4. Time scales t a and t' a and their ratio vs a = T/T c — 1 
in the CFS. Data for the correlation scale c CT are calculated 
exactly by using Eq. (131 with A ~ 0.815. Straight lines are 
MCT predictions. 


obtained from the critical exponents estimated above are 
Ai ~ 1.58, B\ ~ 0.812 and B ~ 1.32. We shall use these 
values to test in a self-consistent manner the MCT pre¬ 
dictions: any observed deviation would mean that MCT 
does not hold for CFS either because the A estimated for 
our system is not connected to the critical decay expo¬ 
nents a and b by Eqs. Q and (|8|, or because Eqs. ([5|, 
([6]) and ([7]) just do not apply to CFS. We have explored a 
relatively wide range of temperatures above T c . Figure [3] 
shows the behavior of 4>(f) for some values of T along 
with the MCT predictions for the crossover through the 
plateau. The full line above f c represents the function 
fc + c a g a (t/t*) while the plateau-below function corre¬ 
sponds to f c + c a gb{t/t*). To appreciate the effect of 
leading order corrections we also add the critical decay 
laws (which are obtained by setting A 1 = Bi = 0) as 
dashed lines. We see that, in agreement with the nature 
of MCT corrections, the quality of comparison increases 
as T decreases. In particular, up to the 80% of the corre¬ 
lator shape is accurately reproduced, over a time window 
ranging from two to more than three decades; see Fig. [5] 
No significant improvement is achieved with the exact 
numerical solution of Eq. (J3| [22]. It is interesting to ob¬ 
serve that, even though MCT was not obviously devised 
for describing CFS, the MCT predictions perform even 
better than a recent truncation scheme for facilitated 
master equations; see Fig. 4 in Ref. [23] for a compar¬ 
ison. Thus, it seems that MCT fully captures in a very 
general way the slow dynamics of large scale cooperative 
rearrangements occuring near the glass singularity. 

What makes the observed agreement pretty remark¬ 
able, however, is the absence of any fitting parameter 
in our procedure. In fact, one should keep in mind that 
the prefactors of leading order corrections, A\ and B 1 /B, 
are correlated and have no temperature dependence, i.e., 
they are strongly constrained. Thus, to check the sta- 
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FIG. 5. Persistence 4>(t) vs time t at temperature T in the CFS. Full lines are the asymptotic leading corrections to the critical 
scaling laws predicted by MCT. The line above the plateau at f c = 0.67309 represents the function f c + c a g a (t/t*) while the 
blue line below the plateau corresponds to f c + c a gb{t/t*). For comparison, pure critical decay laws are shown as dashed lines. 


bility of our comparison, we changed the exponent a by 
an amount of ±0.01, that is essentially the uncertainty 
we have on this exponent. We then found that relaxation 
behavior is still well reproduced if all other quantities en¬ 
tering Eq. |5]) are changed according to MCT relations. 
This robustness seems to suggest that prefactors of lead¬ 
ing order corrections are rather optimal from the purely 
data fitting point of view. Nevertheless, we notice that 
very close to T c (below T = 0.51), the theoretical value 
of c a need to be increased by a small factor (up to 5%) in 
order to accurately describe the intermediate stage of a 
relaxation. This discrepancy is presumably due to criti¬ 
cal finite-size corrections which are hardly accounted for 
in a MCT description. The exact solution of a related 
cooperative facilitated system should help to settle this 
issue [22]. Since static correlations are absent here, we 
expect that a suitable MCT formulation of facilitation 
dynamics is related to those studied in Refs. |25 ] [26 ] , 

Conclusions. - We have shown that the relaxation 
behavior of cooperative facilitated spin models is accu¬ 
rately described, in a relatively wide range of tempera¬ 
tures, by the MCT predictions for the asymptotic leading 
corrections to critical decay laws. These results support 
the idea of glassy universality and put on a solid quanti¬ 
tative ground the correspondence between MCT and CFS. 
At the same time new challenges arise. 

It is well known that the dynamics of disordered p-spin- 
like systems is also exactly described by MCT which is 
the key to the celebrated random-first order transition 
(RFOT) scenario for the glass transition [27; 128]. RFOT 
is based on a nontrivial phase structure characterized 
by a one-step replica symmetry breaking and a possi¬ 
ble fractal free-energy landscape lying deeper in the glass 
phase [25]. In the CFS, instead, thermodynamics plays no 
role: the configuration space breaks into several compo¬ 
nents that cannot be connected by a sequence of allowed 


spin flips. The duality between kinetic and static rep¬ 
resentations seems therefore irreconcilable. This would 
imply, for example, that the dynamical inverse problem 
of inferring the interactions among system units starting 
solely from time correlation data is undecidable. How¬ 
ever, in view of the double connection (between MCT and 
CFS on one hand, and between MCT and RFOT on the 
other) one is led to conclude that a map connecting CFS 
and RFOT must necessarily exist [50] , Extending the 
mapping below the Kauzmann temperature would obvi¬ 
ously require extra thermodynamic ingredients which are 
arguably provided by suitable energetic interactions, as 
shown in Ref. |31j . Thus, the existence of three a priori 
distinct frameworks that give the same dynamical scaling 
laws, while certainly surprising, should not be considered 
as a mere coincidence, but rather as a genuine signature 
of the universality of glassy relaxation captured by these 
frameworks. One advantage of the CFS is that the glass 
formation can be interpreted in real space as the forma¬ 
tion of a bootstrap percolation backbone of permanently 
frozen spins. Much less clear is how to extract the val¬ 
ues of the critical exponents a and b from this geomet¬ 
rical structure, since no spin is permanently frozen on 
the liquid side of the glass transition. So, at the moment 
this important issue remains elusive in this framework as 
well as in other ones. Another perspective suggested by 
the CFS is that the ergodicity-restoring activated events 
in finite dimensions should be interpreted as the ana¬ 
log of the finite-volume metastability effects which are 
known to transform the bootstrap percolation transition 
in a (sharp) crossover [32]. Finally, it would be crucial 
to extend the above investigation to finite dimensional 
facilitated systems having a glass transition (see, e.g., 
Ref. [33]). The exploration of these problems is left to 
future works. 
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